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A  pair  of  masses  or  opposite-sign  charges  re¬ 
leased  from  rest  will  move  directly  toward 
each  other  under  the  action  of  the  inverse- 
distance-squared  force  of  attraction  between  them.  An 
exact  expression  for  the  separation  distance  as  a  func¬ 
tion  of  time  can  only  be  found  by  numerically  invert¬ 
ing  the  solution  of  a  differential  equation.  A  simpler, 
approximate  formula  can  be  obtained  by  combining 
dimensional  analysis,  Kepler’s  third  law,  and  the  fa¬ 
miliar  quadratic  dependence  of  distance  on  time  for 
a  mass  falling  near  Earth’s  surface.  These  exact  and 
approximate  results  are  applied  to  several  interesting 
examples:  the  flight  time  and  maximum  altitude  at¬ 
tained  by  an  object  fired  straight  upward  from  Earth’s 
surface;  the  time  required  for  an  asteroid  of  known 
starting  position  and  speed  to  cross  Earth’s  orbit  if  it  is 
bearing  toward  the  Sun;  and  the  collision  time  of  two 
oppositely  charged  particles  starting  from  rest. 

Problem  Statement 

Suppose  that  the  first  particle,  whose  motion  is  of 
primary  interest,  has  mass  m  and  charge  q,  while  the 
second  particle  has  mass  M and  charge  -Q.  To  ensure 
that  the  electrical  force  between  them  is  not  repulsive, 
assume  that  Qq  ^  0.  (Note,  however,  that  q  can  be 
positive,  negative,  or  zero.  In  the  last  case,  the  force  of 
attraction  is  solely  gravitational.) 

We  will  analyze  the  problem  in  the  center-of-mass 
frame  of  reference  and  assume  that  the  two  particles 
constitute  a  bound,  isolated  system.  In  other  words, 
the  total  mechanical  energy  is  negative  relative  to  the 
usual  potential  reference  at  infinite  separation,  and 


the  particles  always  have  equal  and  opposite  linear 
momenta.  Consequently,  there  is  some  maximum 
distance  Tq  between  the  two  particles  at  which  point 
they  are  both  instantaneously  at  rest.  Define  the  time 
of  that  event  as  t  =  0.  It  corresponds  to  the  initial  sepa¬ 
ration  between  the  particles  if  they  start  from  rest;  it 
represents  their  turning  points  if  the  particles  start  out 
moving  away  from  one  another;  and  it  represents  a 
point  in  the  extrapolated  past  if  they  are  initially  mov¬ 
ing  toward  each  other. 

To  summarize,  without  loss  of  generality  we  sup¬ 
pose  that  two  particles  1  and  2  are  released  from  rest 
at  time  t=0  when  they  are  a  center-to-center  distance 
rg  apart.  Define  the  coordinate  origin  to  be  at  the  sta¬ 
tionary  center  of  mass  of  the  system.  The  two  particles 
attract  and  begin  to  move  toward  each  other.  Their 
radial  separation  and  relative  speed  at  any  later  time 
f  (before  they  collide  with  each  other  at  time  at  the 
origin  if  they  are  of  point  size)  are  r  and  ii,  as  described 
in  Fig.  1 . 

The  attractive  force  between  the  two  particles  is  the 
sum  of  the  gravitational  and  electrostatic  forces,  with 
magnitude  A  =  y  /r^,  where  y  =  kQq  +  GMm.  (Fiere  k 
is  the  Coulomb  constant  and  G  is  the  universal  gravi¬ 
tational  constant.)  In  practice,  one  of  the  two  forces 
will  usually  dominate  and  we  can  neglect  the  other, 
but  it  is  just  as  easy  to  be  general  and  include  both. 
Newton’s  second  law  describing  the  magnitude  of  the 
relative  acceleration  a  =  of  the  particles  toward 

each  other  is' 
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Fig.  1.  Particle  1  of  mass 
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Fig.  2.  A  highly  elliptical  two-dimensional  orbit  of  a 
planet  of  mass  m  shown  at  the  initial  instant  when  it 
is  located  at  aphelion  with  respect  to  a  sun  of  mass  M. 


where  the  reduced  mass  of  the  system  is  fj,  = 
mM  I {m  +  At),  which  simplifies  to  fj,  ~  m  for  the 
common  case  of  Ad  »  m.  The  signs  of  the  kinematic 
variables  have  been  chosen  so  that  r,  tj,  and  a  are  all 
positive  for  0  <  t  <  Consequently  tj  =  —dridt  since 
r  decreases  with  t  during  this  time  interval,  while  a  = 
dvidt  =  -cP-rldf-  because  tj  increases. 

We  can  now  succinctly  state  the  problem;  Find  the 
distance  r  between  the  particles  as  a  function  of  time 
t.  An  exact  solution  for  the  inverse  relationship  r(r) 
has  been  published  in  several  places^"^  and  is  reviewed 
in  the  appendix.  The  purpose  of  the  present  paper  is 
to  develop  an  approximate  formula  using  ideas  that 
would  be  more  appropriate  in  an  introductory  course. 
Not  only  is  the  resulting  expression  for  r(r)  simpler, 
but  it  is  analytically  invertible  to  give  r{t),  in  contrast 
to  the  exact  expression  that  is  not. 

Three  steps  are  used  to  develop  the  approximate 
formula.  First,  dimensional  analysis  quickly  gives  a 
good  estimate  for  the  collision  time  t^.  of  two  point 
particles  that  begin  at  rest.  Second,  Kepler’s  third 
law  determines  the  missing  numerical  prefactor  that 
dimensional  analysis  cannot  give.  Third,  the  familiar 


kinematic  result  for  an  object  falling  in  Earth’s  surface 
gravity  motivates  the  final  expression.  Instructors 
need  not  follow  all  three  steps.  For  example,  one  could 
skip  directly  to  the  second  step  to  find  the  exact  col¬ 
lision  time.  Alternatively,  one  could  omit  the  second 
step  and  simply  state  without  proof  the  small  numeri¬ 
cal  correction  to  the  formula  for  the  collision  time 
found  in  the  first  step.  A  third  possibility  would  be 
to  skip  the  third  step  and  restrict  consideration  to  the 
one-dimensional  collision  of  two  point  particles  start¬ 
ing  from  rest. 


Step  One — Approximate  Solution  for 
the  Collision  Time 


Equation  (1)  implicitly  relates  the  initial  separation 
Tq  to  the  collision  time  t^.  in  terms  of  the  quantity 
Y  //u  that  has  SI  units  of  rrd’ls^.  To  convert  this  quan¬ 
tity  into  one  with  units  of  time,  we  must  divide  it 
by  the  cube  of  and  take  the  reciprocal  square  root. 
Therefore,  we  immediately  predict  that 


(2) 


assuming  both  particles  have  negligible  size  relative 
to  their  initial  separation.  (In  contrast,  an  example 
problem  is  solved  later  in  this  article  where  the  second 
particle  is  taken  to  be  the  Earth  and  its  radius  is  ac¬ 
counted  for.)  As  shown  below,  this  estimate  for  is 
only  10%  smaller  than  the  exact  answer  and  so  it  is 
an  effective  illustration  of  the  power  of  dimensional 
reasoning.^  For  the  case  of  two  uncharged  particles 
(so  that  7  =  GAtni),  when  the  second  particle  is  much 
more  massive  than  the  first  (so  that  [i  —  m),  then  Eq. 
(2)  becomes 


(3) 


Step  Two — Numerical  Prefactor  for 
the  Collision  Time 


The  elliptical  orbit  of  a  planet  (of  mass  rri)  about 
a  sun  (of  mass  Ad  »  rri)  with  semi-major  axis  a  has  a 
period  squared  of 


T 


GM 


(4) 


according  to  Kepler’s  third  law,  as  can  be  derived  at 
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an  introductory  level  using  conservation  laws  and 
elementary  properties  of  ellipses.^  Now  consider  an 
elliptical  orbit  in  the  limit  that  the  eccentricity  e—*\. 
Then  the  sun  (located  at  one  of  the  elliptical  foci) 
moves  to  an  end  of  the  highly  elongated  orbit,  as  illus¬ 
trated  in  a  letter  to  the  editor^  by  a  sequence  of  orbits 
for  which  e  is  increased  from  0  to  1 .  To  make  contact 
with  the  present  paper,  let  us  suppose  that  the  planet 
starts  at  aphelion,  where  its  radial  velocity  is  zero.  In 
that  case,  as  Fig.  2  makes  clear,  the  planet  and  sun  are 
initially  separated  by  the  major  axis  so  that  Vq  =  la.  To 
complete  one  orbit,  the  planet  first  has  to  move  to  the 
sun’s  position,  which  by  definition  takes  the  collision 
time  t^,  and  then  return  to  aphelion  in  the  same  time. 
Thus,  the  orbital  period  is  r  =  2t^.  Substituting  this 
result  and  a  =  r^H  into  Eq.  (4)  gives 


t 


c 


K  I 

2\2GM' 


(5) 


Consequently  the  required  numerical  prefactor  in 
front  of  Eq.  (3)  is  2~^^^  vr  ~  1.11.  Eikewise  correct¬ 
ing  the  more  general  expression  in  Eq.  (2)  leads  to 
the  exact  result  for  two  point  particles. 


2]j  2/' 


(6) 


Noting  that  the  relative  separation  between  the  two 
particles  decreases  by  rg  in  a  time  of  it  proves  conve¬ 
nient  in  the  subsequent  analysis  to  recast  the  position 
and  time  variables  in  dimensionless  form.  Define  the 
normalized  radial  separation  as 


r 

R^-  ij) 

'b 

and  the  normalized  time  as 


(8) 

C 

Note  that  =  1  at  the  initial  point  7’=  0,  and  that 
point  particles  would  reach  the  origin  7?  =  0  at  T  =  1 . 

Step  Three — Approximate 
Relationship  Between  the  Separation 
and  Time 

Suppose  that  a  particle  of  mass  m  is  dropped  near 
Earth’s  surface  (so  that  rg  is  approximately  equal  to 
Earth’s  radius  rg)  and  falls  a  small  distance  Ar  «  rg  to 


a  new  radial  position  r=rg  —  Ar  in  a  time  t.  Elemen¬ 
tary  kinematics  tells  us  that 

Ar  =  —gt^,  (9) 

where  g  =  GMIr\  =  9.8  N/kg  is  Earth’s  surface 
gravitational  field  strength  (with  M  the  mass  of  the 
Earth).  Use  Eqs.  (7)  and  (8)  to  re-express  the  fall  dis¬ 
tance  and  time  in  normalized  form  as 


R  =  \  —  X  where 


A  r 

X  = - <<1 


(10) 


and 


where 


2pg 


(11) 


from  Eq.  (5).  Using  these  variables,  Eq.  (9)  can  be  re¬ 
written  as 


x  =  -r^ 


where  n  = 


(4) 

IT 


T.62.  (12) 


Multiplying  this  equation  for  Xhy  n  and  adding  and 
subtracting  unity  gives 


\-{\-nX)  =  T^. 


(13) 


We  want  to  extend  this  result  to  hold  for  values  of 
X  that  approach  1 .  The  expression  in  parentheses 
looks  like  the  first  two  terms  in  the  binomial  expan¬ 
sion  o£{\  -X)”  =  R",  using  Eq.  (10)  to  relate  A  to 
R.  Making  this  identification  in  Eq.  (13)  leads  to  the 
following  hypothesized  expressions  relating  the  nor¬ 
malized  separation  and  time, 

1-7?”  ^  r(7?)« Vi-7?” 


or 


R{T)^ 


(14) 


Note  that  these  expressions  fit  the  two  endpoints  R=\ 
at  T  =  0,  and  7?  =  0  at  T  =  1 . 

Equation  (14)  is  plotted  in  red  in  Eig.  3.  Eor  com¬ 
parison,  the  exact  result  from  Eq.  (30)  in  the  appendix 
is  graphed  in  blue.  The  two  curves  overlap  almost  per¬ 
fectly  across  their  whole  range!  The  agreement  can  be 
slightly  improved  by  rounding  off  n  to  the  value  1.6, 
which  is  more  convenient  to  work  with  than  the  exact 
irrational  value;  in  that  case,  the  mean  discrepancy^*^ 
between  the  exact  and  approximate  values  of  T{R)  is 
about  a  quarter  of  a  percent. 
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Fig.  3.  Approximate  (red)  and  exact  (blue)  plots  of  the 
separation  of  the  two  particles  as  a  function  of  time  in 
dimensionless  form. 


Three  Example  Problems 

1.  Suppose  a  particle  is  launched  upward  from  the 
surface  of  the  Earth  with  a  speed  "Uorbit  equal  to  that 
of  a  low-Earth-orbit  satellite.  What  is  the  maxi¬ 
mum  height  that  it  reaches  above  Earth’s  surface 
and  how  long  does  it  take  to  get  there?  ^  ^ 

The  speed  of  orbit  of  a  satellite  of  mass  m  just  above 
the  surface  of  the  Earth  (ignoring  its  atmosphere  and 
topography)  is  obtained  from  Newton’s  second  law  as 

te  (15) 

^orbit  =  =7.9  km/s  , 

where  M  and  ^e  are  Earth’s  mass  and  radius,  respective¬ 
ly.  Since  the  velocity  on  the  way  up  is  equal  and  op¬ 
posite  to  the  subsequent  velocity  back  down  from  the 
peak  height,  we  can  substitute  r  =  ^e,  D  =  tJorbit>  und 
7  //u  =  GM  =  gr\  into  Eq.  (25)  to  obtain  the  maxi¬ 
mum  radial  position  tg  =  2rE,  i.e.,  the  particle  attains  a 
maximum  height  of  one  Earth  radius  (6370  km)  above 
the  surface.  It  then  falls  back  to  the  ground,  at  which 
point  R=  r^lrQ  =  112.  Substituting  this  value  into  Eq. 
(30),  one  then  finds  that  the  fall  (or  equivalently  the 
rise)  time  is 


t 


=  34.5  min 


(16) 


since  T  =  where  ^  according  to  Eq. 

(5).  For  comparison,  Eq.  (14)  predicts  a  time  of 


(17) 


using  the  rounded-off  value  of  n.  The  numerical  val¬ 
ues  of  Eqs.  (16)  and  (17)  disagree  by  less  than  1  s. 


2.  A  new  asteroid  is  discovered  that  is  heading 
straight  for  the  Sun.  It  is  initially  observed  at  a 
distance  from  the  Sun  of  3  AU  with  a  speed  of  half 
of  the  speed  it  would  need  to  escape  from  the 
Sun  at  that  point  in  its  orbit.  How  much  time  will 
elapse  until  the  asteroid  crosses  Earth’s  orbit? 

Eet  f;  =  4.5  X  10'^  m  and  =  1.5  X  10^^  m  be  the 
initial  and  final  radial  positions  of  the  asteroid,  respec¬ 
tively,  and  let  rg  denote  the  maximum  distance  from 
the  Sun  that  the  asteroid  would  attain  if  it  were  travel¬ 
ing  directly  away  from  the  Sun  with  an  initial  speed  of 
/2.  Also  let  =  rj Tq  and  R{=  rflr^  .  The  initial 
speed  of  the  asteroid  (of  mass  m)  is  obtained  from  en¬ 
ergy  conservation  as 


1 

—  mv 
2 


2 

esc 


GMm 

- =  0  ^ 

4 


H=-^^esc=. 


\GM 
'  2r: 


=  12.1  km/s , 


(18) 


where  M  is  the  mass  of  the  Sun.  Next  we  use  Eq. 
(25)  with  7  //u  =  GM  to  find  the  reciprocal  of  the 
asteroid’s  radial  turning  point. 


1  _  1  vf  _  3 
Tq  T;  2GM  4r; 


(19) 


with  the  help  of  Eq.  (18).  Consequently,  =  3/4 
and  R^=  1/4.  Thus  Eq.  (30)  gives 


Ar 


T -h  = 


-1  1  -1 
cos  . - cos 

u 


=  — J - =  173  days. 


using  Eqs.  (5)  and  (8).  In  comparison,  Eq.  (14) 
implies  that 


Ar  =  27r 


.^^Ul-0.25^-^  -Vl-0.75 
^GM 


1.6 


(21) 
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(24) 


and  the  quantity  in  square  brackets  is  nearly  equal 
to  1/3,  in  agreement  with  the  result  of  Eq.  (20).  If 
such  an  asteroid  really  were  discovered  tomorrow,  we 
would  therefore  have  about  half  a  year  to  prepare  if 
the  Moon  or  Earth  were  to  be  impacted. 

3.  An  electron  and  a  positron  start  out  at  rest  a  dis¬ 
tance  fg  =  10  nm  apart.  How  much  time  will  it  take 
for  them  to  annihilate  each  other? 

The  two  particles  have  equal  and  opposite  charges 
of  magnitude  q  =  Q  =  e,  where  e=  1.6  X  10“^^  C. 
Furthermore  they  have  equal  masses  m  =  M  =  9.11 
X  10“^^  kg  so  that  the  reduced  mass  is  ^  =  ml2  .  The 
electrostatic  force  completely  overwhelms  the  gravita¬ 
tional  force  between  them  and  thus  7  =  ke^.  Equation 
(6)  now  becomes 

•.■iff-’" 

They  disappear  in  less  than  a  blink  of  an  eye! 


1  2  7  n  T 

—  /j-v - =  0 - , 

2  r  Tq 

which  can  be  rearranged  as 


V  = 


^---1 

i  R 

(25) 


The  square  root  is  always  real  because  r  ^  rg.  Note 
that  the  speed  becomes  infinite  as  r  — ^  0,  but  in  the 
real  world  that  would  be  prevented  by  the  finite  sizes 
of  the  two  particles. 


II.  Relationship  Between  the  Separation  and 
Time 

Substituting  v  =  —dridt  into  Eq.  (25),  the  variables 
r  and  t  can  be  separated  and  rewritten  in  dimensionless 
form  using  Eqs.  (6)  through  (8)  to  get 


dR 


(26) 


Appendix:  Exact  Relationship 
Between  the  Separation  and  Time 

Substituting  a  =  -d^rldR  into  Eq.  (1)  results  in  a 
second-order  differential  equation  for  r{t)  that  can  be 
solved  in  two  steps.  First  we  find  the  relative  speed  v 
of  the  particles  as  a  function  of  their  separation  r ;  then 
we  integrate  that  to  determine  the  time  T  required  for 
the  particles  to  approach  each  other  to  within  some 
distance  R  in  dimensionless  form. 


I.  Relative  Speed  as  a  Function  of  the 
Separation 

In  the  center-of-mass  reference  frame,  the  two 
particles  have  equal  and  opposite  linear  momenta  of 
magnitude^  mV\  =  Mxs2  =  [iv  =  p,  so  that  their  total 
kinetic  energy  is 


2^,2  ( 


P  H 

^  P  - 

fJL  V 

2m 

2M 

2 

1  1 

- 1 - 

m  M 


1  2 

=  —  av  . 

2 


(23) 


The  potential  energy  corresponding  to  the  force  of  in¬ 
teraction  in  Eq.  (1)  is  -7/r .  Conservation  of  mechani¬ 
cal  energy  for  the  system  of  two  particles  between  the 
initial  and  an  arbitrary  final  state  therefore  becomes 


The  lower  limits  correspond  to  the  initial  condition 
=  1  at  7’=  0.  The  time  integral  on  the  left-hand 
side  is  equal  to  T.  One  can  solve  the  radial  inte¬ 
gral  by  making  a  change  of  variable.  By  trial  and 
error,  one  is  led  to  the  trigonometric  substitution 
R  =  cos^  Q,  whose  virtue  is  that  it  lends  itself  to  the 
Pythagorean  identity 

sec^  6-1  =  tan^  6.  (27) 

Substituting  this  identity  and  the  differential  dR  = 

—2  cos  6  sin  ddd  into  Eq.  (26)  leads  to 

e  0 

— T  =  ^ 2cos^  9 do  =  y  (1  +  cos 20) <5/0,  ^28) 

0  0 

using  the  cosine  double-angle  formula  and  noting 
that  0  =  0  when  7?  =  1 .  Integrating  gives 

IT  ^  „  sin  20  ^  ^  ^ 

—T  =  9-\ - ^ —  =  0  +  cos0sm0  (29) 

=  0  +  cos  O^Jl  —  cos^  0 

with  the  help  of  the  sine  double- angle  formula  and 
another  Pythagorean  identity.  Returning  to  the  origi- 
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nal  variable  R  gives  the  final  result 


-T  =  cos~^  ^  +  jRil- R).  (30) 

2 

This  equation  correctly  predicts  that  T  =  0  at  the  ini¬ 
tial  point  ^  =  1,  and  that  T  =  1  at  7?  =  0.  We  cannot 
analytically  invert  it  to  find  RiT),  but  by  computing 
a  table  of  values  (as  plotted  in  Fig.  3)  we  can  find 
any  desired  result  numerically,  which  points  out  the 
advantage  of  working  with  a  “universal”  (i.e.,  dimen¬ 
sionally  normalized)  form  of  the  equation. 

One  can  show  that  Eq.  (30)  becomes  Eq.  (12)  to 
lowest  order  inXif  Eq.  (10)  holds.  This  result  follows 
by  proving  that  each  of  the  two  terms  on  the  right- 
hand  side  of  Eq.  (30)  is  approximately  equal  to  ^fx  if 
R=  \  —X,  where  2f<<  1.  For  the  first  term,  one  dem¬ 
onstrates  the  inverse  result  yj\  —  X  ~  cos  'Jx  by 
Taylor  expanding  both  sides  up  to  linear  powers  in  X. 

Note  added  in  proof:  S.K.  Foong  [“From  Moon-fall  to 
motions  under  inverse  square  laws,”  Eur.  J.  Phys.  29, 
987-1003  (Sep.  2008)]  has  obtained  other  approximate 
expressions  for  R{T). 
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